Modeling habitat suitability of Dorema ammoniacum D Don. in the rangelands of central Iran

The purpose of this study was to evaluate the predictive accuracy of habitat suitability models, identifying the potential distribution range of Dorema ammoniacum, and its habitat requirements in the rangelands of Yazd province, central Iran. Bafgh, Mehriz and Nadoushan, were three habitats that were identified, and sampling was conducted in each habitat using a random-systematic method. A set of 10 plots were established (at equal distances) along 350 m long 18 transects. Soil samples (two depths: 0–30 and 30–60 cm from 36 profiles) were collected and measured in the laboratory. Elevation, slope, and aspect maps were derived, and climate information was collected from nearby meteorological stations. The habitat prediction of the species was modeled using Logistic Regression (LR), Maximum Entropy (MaxEnt), and Artificial Neural Network (ANN). The Kappa coefficient and the area under the curve (AUC) were calculated to assess the accuracy of the forecasted maps. The LR model for habitat prediction of the studied species in Mehriz (K = 0.67) and Nadoushan (K = 0.56) habitats were identified as good. The MaxEnt model predicted the habitat distribution for the selected species in Bafgh and Mehriz habitats as excellent (K = 0.89, AUC = 0.76, K = 0.89, AUC = 0.98), and in the Nadoushan habitat as very good (K = 0.78, AUC = 0.85). However, the ANN model predicted Bafgh and Nadoushan habitats as excellent and Mehriz habitat as very good (K = 0.87, K = 0.90, and K = 0.63, respectively). In general, in order to protect species D. ammoniacum, the development of its habitats in other areas of Yazd province and the habitats under study in conservation programs should be given priority.


Study area
D. ammoniacum has three major habitats in Yazd Province, Iran, and is distributed in Nadoushan, Bafgh and Mehriz regions (Fig. 1).These maps of studied areas were generated by the researchers using ArcGIS version 10.1 39 .Furthermore, the key specifications of the study regions are elucidated in the subsequent sections.(i) Gazestan Bafgh is in the southwest of Bafgh county (55° 42' to 56° 2' E and 31° 26' to 31° 47' N).The altitude range is 1800-2600 masl, and the average annual precipitation is 150 mm per year.The average annual temperature is 15.5 °C.The soil texture dominating this area is sandy and loamy.The composition of the livestock herd consists of approximately 70% goats and 30% sheep.(ii) Aliabad Ghelgazi Mehriz is in the west of Mehriz county with a geographical range of 54° 5' to 54° 34' E and 31° 10' to 31° 24' N. The altitude range is 1800-2800 masl.The average annual precipitation is 140 mm per year, and the average annual temperature is 13.5 °C.The composition of the livestock herd consists of approximately 70% goats and 30% sheep.(iii) Nadoushan is in the west of Meybod county with a geographical range of 53° 1' to 53° 45' E and 31° 28' to 31° 38' N.This area is primarily plain with an altitude range of 1370-2900 m above sea level (masl).According to the 22-year statistics of meteorological stations in the study areas, the average annual precipitation is 110 mm, and the average annual temperature is 18.5 °C.Soil texture is sandy and sandy-loamy.According to field observations and surveys conducted with local

Data collection
Habitat identification within each region involved a comprehensive approach, combining landform and geological maps, field surveys, and observational data to delineate sampling units 30,[41][42][43] .Following this, areas devoid of the target species were randomly chosen near existing presence areas, ensuring similarity in soil composition and ecological conditions.Subsequently, non-presence sites for the studied plant species were selected through a random process.In each of the identified habitats, sampling was performed using a systematic-random method along 18 sampling transects of 350 m.On each transect, 10 plots of 2 m 2 were established.The length of sampling transects, and quadrat size was determined based on the density and type of vegetation identified in previous studies 30,[41][42][43] .Two soil profile samples were collected along each transect.The depth of the profiles was determined according to the effective rooting depth of the studied species.As the maximum rooting depth was identified as 60 cm on average 30 , two soil samples were collected from 0-30 cm as the surface soil and 30-60 cm as the subsurface soil.A total of 36 soil samples were collected from each habitat.Then, several soil physical and chemical properties were measured in the laboratory, including texture (Bouyoucos hydrometer method), saturation moisture (SP), electrical conductivity (EC) (saturated extraction method), acidic and alkaline (pH) (potentiometric method), organic carbon (OC) and organic matter (OM) Walkely and Black 44 , soluble solutes including Na and K (flame photometer method), Ca, Mg, K, Cl-, CO 3 -and HCO 3 -(titration method) and SO 4 -(spectrophotometer measured in the laboratory) 30,32 .The geographical location of the sampling points was recorded using a hand-held Global Positioning System (GPS).Soil factor maps affecting the habitat of D. ammoniacum were generated through Kriging interpolation using ArcGIS version 10.1 39 .Soil factor coefficients were derived from regression analysis and then integrated into the environmental layers.The desired parameters for each sampling point were then extracted from slope, aspect, and altitude maps by the researchers using Arc-GIS version 10.1 39 .Meteorological data (average annual precipitation and temperature) were collected from the meteorological stations from the study region.Total data were organized in Excel according to different habitats, followed by a check for normality using Kolmogorov-Smirnov and Anderson Darling tests.If normalization was necessary, methods such as log, square-root, and Box-Cox transformation were applied; otherwise, nonparametric geostatistical methods like inverse distance weighting (IDW) were utilized.Prior to model application, www.nature.com/scientificreports/initial multicollinearity among independent factors was assessed, and factors with VIF > 10 were eliminated 45 .Environmental factor maps were generated using ArcGIS 10.1 software, and Table 2 provides details on the total number of variables, their determination method, and corresponding references.

Data analysis
In this study, three machine learning algorithms were employed to predict plant habitat: Logistic Regression (LR), Maximum Entropy (MaxEnt), and Artificial Neural Networks (ANN).These methods offer diverse approaches to modeling and analyzing habitat distribution, each with its own strengths and applications.

Logistic regression (LR) model
Logistic Regression (LR) is a highly valuable method for exploring the relationship between independent variables and a binary response variable, such as the presence or absence of a specific plant species.As a specialized form of multiple regression, LR is particularly adept at analyzing discrete dependent variables.Presence-absence models, like logistic regression, frequently yield ecologically sound relationships, providing valuable insights into species distribution dynamics 61,62 .In this study, the logistic regression model, as represented by Eq. 1, was employed to predict the presence and absence of plant species 63 .In LR model, a dataset comprising 180 points is utilized for modeling.This dataset consists of 90 points indicating species presence and 90 points indicating species absence.Here, plant species served as the dependent variables, while environmental factors acted as the independent variables.Coefficients were assigned to each layer within ArcGIS, facilitating the generation of prediction maps for the studied species across different habitats.These maps provided a continuous probability range from 0 to 1 based on the presence or absence of D. ammoniacum species.The logistic regression model was implemented using SPSS Ver.18 software.
where Y is the probability of occurrence of the species, x are the predictive variables and b are the regression model coefficients.

MaxEnt model
The Maxent model, a species distribution model derived from machine learning, is utilized to forecast the potential distribution of species 64,65 .This model, which estimates species distribution based on presence-only data, has demonstrated practical effectiveness 21,32 and can aid in predicting the detectability of a species [66][67][68] .Incorporating both continuous and categorical data as well as interactions between variables, this method has been shown to create satisfactory results even with limited sample sizes.By determining the optimal distribution function between species presence points and environmental variables using the maximum entropy principle, the MaxEnt model then extrapolates this function across the study area to generate a distribution map of the target species 21,69 .
In this method, first, the environmental layers are evaluated as the location of training data, and then the probability of occurrence of the studied species in the whole region is selected 32,68 .Data of 90 occurrence points from sampling sites were used for model building.In regards to spatial thinning, it's important to mention that the selection of sampling sites and plant presence locations has been accurately designed to prevent false replicates and ensure representative soil samples of suitable spatial distributions based on environmental diversity.The rationale behind the research design was to carefully select sampling sites and locations where the plant species is present in order to avoid redundant data points and ensure a diverse representation of environmental factors.This approach was intended to facilitate a spatially balanced distribution of samples of the target species and soil samples, aligning with the research objectives.
Maps of environmental factors built into ASCII template, and MaxEnt 3.3 software applied for modeling the species distribution.In the algorithm, the number of iterations was set to 5000, with a convergence threshold of 0.00001.A random test percentage of 25 was allocated for test data, and the "Random seed" option was checked for replicate run types.Additionally, a subset was utilized for most of the implemented models, with 10,000 background points, a regularization multiplier of 1, and 15 runs in the replicates field.The output grid format was logistic, and the algorithm factors were set to "auto features", with default values chosen for other options 70 .In the process, 25% of the data were separated for model testing (accuracy assessment) and the remaining data were used for training 67 .The receiver operating characteristic curve (ROC) was obtained to test the efficiency of the MaxEnt model.Additionally, the area under the curve (AUC) value was computed for assessing the precision of the model.AUC presents the overall accuracy of the model 68 .The value range for AUC is from 0 to 1 71 .AUC scores above 0.5 indicate the model performs better than random and above 0.75 are considered acceptable for predictive accuracy 72 , with > 0.8 considered good and 0.9 considered excellent 73 .The AUC statistic presented the model's superiority in recognition among existence and non-existence areas.The statistical value near unity showed a better consistency of the prediction model with the recorded sampling data and reality 74 .We applied the Jackknife test to distinguish the most effective environmental factors 18 .Jackknife of regularized training gain is very important in determining the most important environmental factors affecting species distribution.If this model is implemented separately, it compares the training gain of each variable and then compares it with all environmental variables 10,32 .The probability of species presence ranging from 0 to 1 was estimated 75 .

Artificial neural networks (ANN) model
ANN is a species distribution model that predicts the distribution of a species using presence and absence data.ANN is far from the drawbacks of traditional systems and has the ability to learn and generalize using the examples provided during the training phase.This model is a simulation method and one of the most common classifications for pattern resetting 76 .The advantage of the neural network method is to learn directly from the data, without the need to estimate their statistical characteristics 77 .The neural network is able to find the relationship between a set of inputs and outputs to predict any output corresponding to the desired input and is able to estimate any type of nonlinear function 78,79 .
(1)  80 .These networks consist of input, hidden, and output layers, with the number of neurons in each layer determined by the specific problem being addressed.The configuration of hidden layers and their neurons is typically optimized through trial and error to minimize errors 81,82 .The process of utilizing MLP involves creating, training, and validating the network.During network creation, inputs and outputs are defined, and model performance is assessed using a quadratic polynomial equation (Eq.( 2)) 83 .The training phase involves iteratively adjusting network parameters to minimize the root mean square error (RMSE) using a subset of the data.where x i and x j are any pairs of the independent variables, while z ij is the predicted output from the first layer of neurons in the GMDH algorithm.The described procedure is repeated for a given set of 'n' observation of the m independent variables (x 1 , x 2 , > … > , xm).As the number of input variables is m, the total number of combination of each pair of variables x 1 and x j is m (m − 1)/2.
The processing and creation of spatial change maps resulting from species distribution modeling were conducted by the researchers using ArcGIS version 10.1 39 , employing diverse functionalities such as spatial and 3D analysis, and subsequently presented as spatial maps.
To validate the maps obtained from the model prediction, the presence data that were not analyzed were used as a reference image.Also, in this study, Kappa statistics were used to assess the performance of the models by comparing the probability forecast map with the actual map.To compare the observed and predicted maps, a cut-off point is required to convert continuous probabilities to binary probabilities (presence and absence).In the present study, the presence points predicted in the prediction maps and the presence points in the actual maps were calculated using SPSS Ver.18 software.After determining the optimal threshold, the presence and absence maps of the studied species were classified based on this threshold.Then, the Kappa coefficient was used to evaluate the accommodation between the forecast map and the real map (Eq.( 3)).The accuracy of the forecast maps and real maps were tested in TerrSet 18.31 software 15,32 .

Ethics approval and sampling permissions
All experimental procedures underwent approval from the Review Board of the Department of Range and Watershed Management, Faculty of Agriculture and Natural Resources, University of Mohaghegh Ardabili, Iran.Adherence to applicable guidelines and regulations was strictly observed throughout the study.Notably, no sampling of D. ammoniacum occurred, and there was no harvesting of any part of the plant.Permission for research activities within the study area was granted by the General Department of Natural Resources and Watershed Management of Yazd Province, Iran.Furthermore, the species has previously been identified, and its habitat has been modeled in this research.

Logistic regression
In the Nadoushan habitat, soil parameters emerge as the primary influencers of D. ammoniacum distribution, as indicated by Eq. ( 4).Notably, the presence of D. ammoniacum exhibits a positive correlation with certain soil characteristics, including the percentage of saturated moisture (SP), calcium (Ca), and sodium (Na) within the surface layer (0-30 cm), as well as organic carbon (OC) within the subsurface layer (30-60 cm).Conversely, the species presence shows a negative association with nitrogen (N) levels within the surface layer.According to the specified optimal threshold (0.21), the degree of accommodation of the forecast map with the real map is at a good level (0.556) (Table 5).
Soil parameters are also the most important factors in the distribution of this species in the Bafgh habitat (Eq.( 5)).The presence of D. ammoniacum has a positive relationship with electrical conductivity (EC), pH and OC of the surface layer (0-30 cm), and a negative relationship with EC and OC of the subsurface layer (30-60 cm).According to the specified optimal threshold (0.11), the degree of accommodation of the forecast map with the real map is at a medium level (0.5) (Table 5).
Soil parameters are also identified as the most effective factors in the distribution of this species in the Mehriz habitat (Eq.6).The presence of D. ammoniacum has a positive relationship with K, HCO 3 and SO 4 of the surface layer (0-30 cm), and a negative relationship with Silt, Clay and SO 4 of the subsurface layer.According to the specified optimal threshold (0.35), the degree of accommodation of the forecast map with the real map is at a good level (0.67) (Table 5).The modeled presence and absence area of D. ammoniacum is shown in Fig. 3.

Maximum entropy
The distribution map of species D. ammoniacum obtained using the MaxEnt method is presented in Fig. 2. Figure 2 illustrates the area under the receiver operating characteristic (ROC) curve, commonly referred to as AUC.These AUC values facilitate straightforward comparisons of model performance and are instrumental in evaluating various MaxEnt models 15 .The accuracy of the prediction model based on the classification is at an acceptable level for Bafgh habitat (AUC = 0.76) and the Nadoushan habitat (AUC = 0.85).For the Mehriz habitat, the accuracy of prediction is excellent (AUC = 0.98).
In Bafgh, the results of the Jackknife method showed that the factors of OM, OC, and N of the surface layer and EC and sand of the subsurface layer are the most important and effective factors.At Mehriz habitat, Na of subsurface layer, Ca, and Silt of surface layer, and slope have a significant relationship with the distribution of D. ammoniacum.At the Nadoushan habitat, the most important variables with the largest share in the model are: N, OC, OM of surface layer, and N, OM of subsurface layer.The results of the accuracy assessment of the models (Kappa coefficient) are shown in Table 5.The final prediction maps were based on two levels of presence (1) and absence (0) of plant species (Fig. 3).

Artificial neural network
Various parameters of the ANN, including the momentum training module, the axon transfer tangent function, and the number of hidden layer neurons and processing elements, were tested to design an appropriate network.The trial and error method was employed to refine the neural network model 83 .Th e results of the potential distribution of D. ammoniacum in the studied habitats by ANN indicate that this species may be present in wider areas.Table 3 shows the results of the implemented models for each habitat.The model was executed using a sigmoid function with 10,000 replications.Accuracy rate was calculated for Bafgh and Nadoushan habitats as 96.43% and 93.33% and for Mehriz habitat as 100%.The forecast map model was prepared after training (Fig. 3).
Based on the results of the ANN, the most important variables affecting the distribution of species D. ammoniacum in Bafgh habitat are altitude, N of surface layer, and EC, N, OC, pH of subsurface layer.In Mehriz habitat, the most important variables affecting the distribution of the studied species are pH of the subsurface layer, Ca, N, SO 4 of the surface layer, temperature, altitude, precipitation, and slope.In the Nadoushan habitat, the most important variables affecting the distribution of the studied species are Clay, K, N of surface layer, and N of the subsurface layer (Table 4).
In general, the range of changes of the most important variables affecting the distribution of species D. ammoniacum across all three habitats and based on all three models used, are detailed in Table 4.The environmental factors that consistently impact species distribution within each habitat are highlighted in bold within the table.The Maxen, Maxent, and ANN models emerged as the most accurate predictors for the Bafeq, Mehriz, and Nadoushan habitats, respectively, as highlighted in Table 4.
The degree of agreement of the prediction maps obtained from the neural network model with the terrestrial reality in the studied habitats of D. ammoniacum is presented in Table 5.Based on the results obtained from the evaluation of the models with Kappa coefficient, the degree of agreement of the predicted map with the terrestrial reality related to Bafgh habitat was excellent.The rate of adaptation was good for Mehriz habitat and excellent for Nadoushan habitat.
The LR, MaxEnt, and ANN models were used to determine the areas classified as suitable and unsuitable for D. ammoniacum within each habitat, as shown in Table 6.The results of the modeling framework offer valuable  insights into the potential presence of D. ammoniacum in the study area, facilitating the implementation of management strategies for its effective long-term conservation.

Discussion
Endemic and rare species have acquired top preference for conservation because of their extinction risk.We used three methods of LR, MaxEnt, and ANN for conservation of D. ammoniacum, a threatened and endemic species of Yazd province.An assessment of the model efficiency in this study showed that the accuracy of the models is different in predicting the presence or absence of species studied.Although, it should be noted that no modeling method is able to perform best in all conditions and the selection of the optimal modeling method, in addition to statistical considerations, is a function of factors such as the purpose of modeling, the ecological characteristics of the species, the type of data available and the interpretability of the results from an ecological point of view 22 .In general, modeling, simulating, and mapping potential habitats for these species can aid in conservation planning and management.Salam et al. 1 reported similar results.In this study, the modeled maps showed good predictive accuracy (Kappa coefficient = good, very good, and excellent), and the environmental factors were identified as important determinants of the D. ammoniacum distribution.According to the Kappa coefficient, the best and most accurate model in predicting species distribution in the study area was MaxEnt, then ANN  the plant distribution and organic matter was 0.85 84 .Soil OM and OC are the most important factors affecting species diversity and should be considered in managing the promotion of species diversity 85 .Additionally, in this study, soil N was also identified as one of the effective factors on the distribution of the studied species.One of the reasons for this can be that N is a very effective element in the expansion and presence of plant species and has the highest amount of element absorption in plants.In other words, the effectiveness of these elements is due to their high mobility 86 .Also, N, OM, and OC are indicators of soil quality (Jannat Babaei et al. 2019).
Habitat requirements of D. ammoniacum indicated that the highest presence probability of this species occurred in soil with high OM (0.42-1.39%),OC (0.17-0.54%),N (0.021-0.07%) of surface layer and sand (47-84%) of subsurface layer.In Mehriz habitat, the accuracy of modeled maps using MaxEnt method was the highest among the methods studied.The soil and physiographic properties were identified as the most effective factors in the distribution of species D. ammoniacum.Sodium (Na) of subsurface layer (30-60 cm), slope, Ca, and Silt of surface layer (0-30 cm) had the greatest effect in the presence of the species.The probability of the presence of species D. ammoniacum in this habitat is inversely related to the increase in soil salinity or sodium; And with increasing soil salinity, the probability of the presence of the studied species decreases.It can be concluded that saline soil cannot meet the ecological needs of this species.The tolerance range of the studied species is from 0 to 10.9.In soils with high salinity, due to high Na ions, poisoning may occur in the plants.Also, large amounts of NaCl in soil may affect plant growth, germination, and production 16 .As was found in this study Ghorbani et al. 30 , also reported a direct relationship between the predictions of the distribution of Leucanthemum vulgare Lam., with Na amount.In studying relationship between soil parameters and the distribution of species D. ammoniacum and Rheum ribes in rangelands of Baghedar region, Iran.Na was also identified as the most effective parameter 87 .Another influential factor determining the distribution of the studied species is Ca of surface soil layer.In investigating the relationship between soil, and vegetation in saline lands of Qom province (Iran), K and Ca were identified as the most influential factors in the community of Halanthium rariflorum 88 .Also, in order to determine the relationship between soil properties and vegetation in Baghdar rangelands in Bafgh, Iran, Ca had the greatest effect on the separation of D. ammoniacum and Rheum ribes 87 .The results showed that slope is another important factor among other variables.The probability of presence of species D. ammoniacum in this habitat is inversely related to the increase of slope and with increasing slope, the probability of presence of the studied species decreases so that it can be concluded that high and mountainous areas reduce the probability of presence of this species.In confirmation of this research finding, it is reported that slope is the only non-soil factor that has a moderate correlation with species composition 89 .Also, slope percentage is one of the most important factors that affect plant moisture.The effect of slope was also recognizable in the depth of the soil and thus, in plant establishment 15 .The relationship between slope percentage and plant species distribution has been emphasized in other studies 22 .The amount of silt had a significant effect on the occurrence of this species especially D. ammoniacum.In a specific climatic area, soil texture has a greater impact on successful plant growth and regeneration than soil chemical www.nature.com/scientificreports/fertility and is one of the features in the growth and establishment of plant species 68 .Studies by other researchers have emphasized the relationship between the percentage of silt and soil texture in the distribution of plant species 17,23,89 .The habitat requirements of D. ammoniacum indicated that the highest presence probability of this species occurred in soil with high Na (0.82-10.9%),Ca (2.1-9.5%), and Silt (4-18%) content.In Nadoushan habitat, the accuracy of the ANN modeling method was the highest.Based on the results of sensitivity analysis, the most important variables affecting the distribution of the studied species were clay, K and N of surface layer.In this habitat, the percentage of clay at the first soil depth is effective in dictating the distribution of this species.So that the probability of the presence of this species has an inverse relationship with the percentage of soil clay and with increasing soil clay the probability of the presence of the studied species decreases.The soil textures in this habitat are sandy, sandy-clay and sandy-loam.In fact, the relationship between clay and vegetation confirms the relationship with soil texture (percentage of clay, sand, and silt).The impact of soil texture in increasing the habitat suitability of plant species has been confirmed in other studies 14,90 .Another effective factor determining the distribution of the studied species was soil N. Nitrogen is important in soil OM and plays an effective role in soil fertility 18 .The importance of N on plant distribution has been emphasized by other researchersv 15,18,22 .Another influential variable in this study is K.Other researchers have concluded that soil nutrients are a key factor in determining plant distribution 43,90 .Studying the habitat requirements of D. ammoniacum indicated that the highest presence probability of this species occurred in soil with high N (0.007-0.067%), clay (6-18%), K (0.9-14.6%), and lime (12-15%) content.

Conclusions
D. ammoniacum is an endemic and threatened plant species in the rangelands of Yazd province, central Iran.It is facing imminent threat due to over-grazing, fragile habitat, construction of roads and over-exploitation for local use, and requires immediate conservation.In Bafgh and Mehriz habitats, among the modeling methods used, the MaxEnt method based on Kappa coefficient had better predictive performance.The MaxEnt method had a better performance because it estimated the probability of occurrence of the studied species based on the limitations related to environmental factors and by the model of MaxEnt.Accordingly, in Bafgh habitat, the variables of OC, OM and N had the highest share in the occurrence of the studied species and in Mehriz habitat, the variables of Na, Ca and soil Silt percentage had the highest share in the presence of the selected species.Therefore, the studied species prefers soil habitats with higher OM and OC, Na percentage and consequently lower salinity, lower silt content and lighter texture.According to the available information, the northern areas of Bafgh habitat and the southern areas of Mehriz habitat are more suitable for the establishment of the studied species.In Nadoushan habitat, among the modeling methods used, the ANN method based on kappa coefficient had better predictive performance.The ANN method performed better because of its ability to recognize complex nonlinear relationships.According to this method, which has produced the most accurate results, the variables of clay percentage, K and soil N had the largest share in the occurrence of the studied species.According to the available information, the southern areas of Nadoushan habitat are more suitable for species establishment.In general, in the present study, we found that MaxEnt and ANN are the best modeling methods.
Currently, the acquisition of accurate information for protective measures and prioritization poses challenges.Documenting the distribution data of the valuable plant D. ammoniacum stands as a crucial step towards addressing these challenges.Predictive maps derived from modeling efforts can aid rangeland experts and managers in better prioritizing actions, such as identifying suitable regions for D. ammoniacum establishment and restoration, and designating conservation units for its preservation.Future studies should extend beyond assessing suitable cultivation areas to exploring changes in habitat size under different climate change scenarios, thus enabling decision-making based on real conditions.Population viability analysis can further enhance conservation efforts by providing insights for effective species preservation.Given the research priority accorded by the General Department of Natural Resources and Watershed Management of Yazd Province, Iran, to investigate and protect D. ammoniacum habitats, the research encountered no significant limitations due to collaborative efforts with various organizations.

Figure 3 .
Figure 3. Modeled maps of D. ammoniacum species using three methods in three habitats (The processing and creation of spatial change maps resulting from species distribution modeling were conducted by the researchers using ArcGIS version 10.1 39 ).

Table 1 .
Characteristics of meteorological stations and climatic variables in the habitats of the study area.

Table 2 .
Environmental variables in modeling habitat suitability of Dorema ammoniacum D Don.

Table 5 .
Optimal presence threshold and degree of accordance of D. ammoniacum species prediction map with terrestrial reality based on three mode.

Table 6 .
Level of suitable/unsuitable area using LR, MaxEnt and ANN models based on thresholds in modeling habitat suitability of Dorema ammoniacum D Don.